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The design of efficient quantum information processing will rely on optimal nonequilibrium transi¬ 
tions of driven quantum systems. Building on a recently-developed geometric framework for comput¬ 
ing optimal protocols for classical systems driven in finite-time, we construct a general framework for 
optimizing the average information entropy for driven quantum systems. Geodesics on the parame¬ 
ter manifold endowed with a positive semi-definite metric correspond to protocols that minimize the 
average information entropy production in finite-time. We use this framework to explicitly compute 
the optimal entropy production for a simple two-state quantum system coupled to a heat bath of 
bosonic oscillators, which has applications to quantum annealing. 
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I. INTRODUCTION 

Though a unifying set of principles describing all 
known nonequilibrium phenomena remains undiscovered, 
many recent developments have illuminated the thermo¬ 
dynamic behavior of small-scale systems. For instance, 
fluctuation theorems valid far from equilibrium have been 
developed in the classical setting mm as well as in the 
quantum regime EM- 

An area of nonequilibrium thermodynamics of partic¬ 
ular interest concerns the operation of small-scale infor¬ 
mation processing systems. The interplay between infor¬ 
mation as a physical quantity and thermodynamics has 
a rich history 15311531 . 

The physics of information processing is of particular 
relevance considering the rapid development of informa¬ 
tion technology and the inevitable approach to compu¬ 
tational limits imposed by physical law |261 EH- Opti¬ 
mization schemes for small-scale information processing 
occurring in finite time will be needed to develop tech¬ 
nology capable of approaching those limits [55J [33]. 

Current research has uncovered techniques to optimize 
thermodynamic quantities arising in small-scale systems 
designed to store and erase classical information 13UH33] . 
including the derivation of a refined second law [35]. This 
research couples with the progress made on the general 
problem of predicting optimal protocols to drive classical 
systems between stationary states with minimal dissipa¬ 
tion [33H3T] . 

In parallel with classical developments, a greater un¬ 
derstanding of optimal processes in the nonequilibrium 
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quantum regime and the efficiency of quantum engines 
has been achieved [42l - i46l . The success of a recently- 
proposed linear-response framework for optimal driving 
of microscopic classical systems [31] calls out for an ex¬ 
tension to quantum systems. 

In the geometric formulation of |34j . a generalized in¬ 
verse diffusion tensor induces a Riemannian manifold 
structure on the space of parameters, and optimal pro¬ 
tocols trace out geodesics of this inverse diffusion tensor. 
This geometric framework is subsequently developed and 
exploited in J33J [35j j3B]. 

In this paper, we extend this work to provide a geomet¬ 
ric framework for computing control protocols optimiz¬ 
ing the average information entropy production [14] . The 
production of entropy is intimately related to the over¬ 
all performance of thermodynamic devices by in some 
sense quantifying irreversibility and providing a bound 
on the availability of useful work. Entropy production 
also provides a useful tool in the analysis of nonequilib¬ 
rium effects. 

By twice measuring the density operator of a system 
interacting with a thermal reservoir at the beginning and 
end of the protocol, the average information entropy pro¬ 
duction may be defined and shown to satisfy a fluctua¬ 
tion theorem m • While there is still debate about the 
best way to define thermodynamic quantities along quan¬ 
tum trajectories, this approach provides an avenue for 
experimental observation M ■ Fortunately, this formal¬ 
ism holds for open quantum systems driven arbitrarily 
far from equilibrium. 

We begin by constructing a general positive semi- 
definite tensor on the space of control parameters for the 
quantum system interacting with the thermal reservoir. 
We assume the dynamics of the system are described by 
a master equation of Lindblad form, arising from an adi- 
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abatic, rotating-wave approximation in a sense we make 
concrete below [47] . 

With the general tensor in hand, we compute optimal 
protocols for a simple two-state system coupled to a ther¬ 
mal bath of bosonic oscillators. The system Hamiltonian 
may be interpreted as describing a spin-4 particle cou¬ 
pled to a time-dependent magnetic field with components 
in the y and ^-directions. This system has applications in 
quantum annealing [47j [48] . For this simple system, we 
demonstrate the existence of null directions of the met¬ 
ric tensor, which correspond to directions in parameter 
space in which there is no overall change in information 
entropy. We derive an approximate expression for the 
optimal overall entropy production. 


II. THE QUANTUM TENSOR 

Our model consists of two distinct components: the 
system under our control and a large collection of quan¬ 
tum degrees of freedom acting as the thermal reservoir. 
Together, the system and the bath degrees of freedom 
evolve unitarily according to the von Neumann equation 
dtPtot = -\\Htot, Ptot\, where 

Htot = H sys ( t ) + Hb + g A a 0 B a . (1) 

a. 

The operator H sys (t) is the system Hamiltonian acting in 
the Hilbert space of system states while Hb is the Hamil¬ 
tonian for the bath degrees of freedom. The interaction 
term g Y) a A a ® B a consists of a (weak) coupling g and 
Hermitian operators A a and B a operating on the system 
and bath Hilbert spaces respectively. 

We are only interested in the time evolution of the 
density operator of the system, denoted pt , which can be 
obtained from p to t by tracing over the bath degrees of 
freedom. We follow the construction of SZI. which ob¬ 
tains a quantum Markovian master equation governing p t 
via an adiabatic, rotating-wave approximation. Specifi¬ 
cally, the derivation of Eq. Q in [47] utilizes the so-called 
“standard adiabatic approximation” 


where r denotes the total evolution time, A = 
min {ei ( t ) — £o(^)} is defined to be the minimum ground 

te[0,r] 

state energy gap of H sysi and 

h= max \(e a {s)\d s H sys (s)\e b (s)}. (3) 

[0,l];a,6 

Here, s = t/r is a dimensionless measure of time. Though 
it is possible to compute higher-order terms in 1/r [47] . 
we will assume that Eq. 0 adequately approximates the 
time-evolution of the quantum system over finite, but 
sufficiently long time scales. 


Assuming a weak coupling g between the system and 
bath degrees of freedom, we have a master equation in 
Lindblad form: 

dtPt = —-j^[H t ,pt] +g 2 C(p t ), 

where 

£(Pt) =^^'ya0{u)( L u,A t )Pt L i,a( t )- 

ctfi co ' 

( 4 ) 

Here, H t = H sys {t) + g 2 H L s{t) consists of the system 
Hamiltonian H sys (t) = H sys (\(t)) and the Lamb shift 
Hamiltonian, which arises through the coupling of the 
system with the thermal reservoir. 

Assume that the system Hamiltonian has time- 
dependent eigenvalues e a (t) with time-dependent eigen- 
kets |e a (f)). Then the operators L WjQ! (f) are defined by 

L u , a {t) = ^ L aba (t) (5) 

,(i)-e a (i)=w 

with 

Lab, a it) = (e a (t)\A a \e b (t))\e a (t))(e b (t)\. (6) 

Furthermore, 

M (7) 

cx.fi co 

where both ^ a pi^) and S a p{io) are Hermitian and can 
be computed from the spectral-density matrix 1:47]. To 
ease notation we will suppress the time dependence of 
the Lindblad operators L WjQ .. 

It is crucial to note that the time dependence of the 
terms defining C arises only through the time dependence 
of the spectrum and eigenkets of the system Hamiltonian. 
Therefore, the time dependence of the Lindblad operator 
C stems from the control parameters A(t). If time ap¬ 
peared explicitly in the terms defining £, we could not 
interpret the approximation developed in this section as 
giving rise to a semi-definite metric on the space of con¬ 
trol parameters. 

In what follows, we denote the control parameter pro¬ 
tocol by A. We assume the protocol to be sufficiently 
smooth to be twice-differentiable. The framework in the 
classical setting is versatile and can handle situations in 
which jump discontinuities are present [33]. Jump dis¬ 
continuities at the end points of the protocol commonly 
arise in optimal finite-time driving processes of classical 
systems [snsuMsiiiis]- For simplicity we only admit 
twice-differentiable protocols though one could in princi¬ 
ple extend this approach to piecewise-continuous A. 

By definition [141 . the average information entropy is 
given by 

(Hj)a = [ dt Tr{(9tp[- lnp + lnp^]}, (8) 

Jo 
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where 


eq 

Pt = 


e ~pH sys (X{t)) 
Tv{e-PHsys (-M0) j- 


(9) 


is the equilibrium distribution defined by the instanta¬ 
neous control parameters A (t) and (3 is related to the 
thermal bath temperature by /3 = (ksT) 1 where ks is 
Boltzmann’s constant. Using Eq. Q, we see 

(£/)a = J dt Trj i[H t ,p t ] + g 2 C(p t )) x 
(-Inpt + lnp| 9 ) 


= 9 2 J dt Tr|>C(p t ) (-lnp t + lnp® 9 ) |, (10) 

where r represents the duration of the protocol. This 
equation follows from the cyclic property of the trace: 

The trace term involving — In p t and the commuta¬ 
tor vanishes because we may permute — In p t into the 
commutator with p t . As for the p e t q term, it too can 
be permuted into the commutator and [H sys , p eq ] = 0, 
naturally. A little more work goes into showing that 


[■ H LS,Pt q } = 0. Since p e t q = 


„-/3 H ays (\(t)) 


Tr{e -/ 3 J bj' s(x(t)) } 

nition, it is sufficient to show that [Hls, H sys {X{t))] = 0 . 

From Eq. ([T]), we need only establish that 
[L'l ta L Ut p,H ay8 (\(t))] = 0. First, note that 


[-Fcj ,/3 i H S ys\ — E (^a |-^/3 |^a) (^61 ? ^si. 

_e b -e a =u> 

= 'y ( (Ca I (Cfo) |e a ) ( e 6 | J 
eb — £a,=Lj 

= ( e a\Ap\e b ){e b -e a )\e a )(e b \ 


e b —e a =u 
= U)L u ^ 


( 11 ) 


where we have suppressed the time dependence for 
brevity. This further implies that 


system and the bath. The system would then evolve in 
time unitarily and the average entropy would vanish. 

This observation allows us to drastically simplify the 
mathematics since the expressions inside the integral 
need only be kept to 0-th order in g 2 . In other words, 
we may compute the evolution of p t using only the von 
Neumann equation 


Z 

&tPt [dd S y S (t), Pt\ - 


(14) 


We wish to approximate Eq (10 1 when the protocol 
duration r is large in the sense of Eq. ([2]). To achieve 
this end we utilize the so-called derivative truncation 
method i.'?3 35; [36], which assumes a specific form for 
the density operator in terms of the equilibrium system 
density operator and the first order derivative of the pro¬ 
tocol A: 


~ e i _i_ A dA “ 

p,*p t +Spx- W - 


(15) 


The Einstein summation convention is assumed here and 
throughout for the index a and the operator 8p\c is Her- 
mitian and traceless. 


Substituting Eq. (151 into the von Neumann equation 


by defi- Eq. ( 14 ) an <4 ignoring derivative terms of order higher 


than first, we obtain equations for the unknown operators 
Sp\a where a indexes the finite set of control parameters. 

These equations are most conveniently expressed in 
terms of the operator basis \uj a (t))(uj b (t)\ where | uj a (t)} 
is an eigenket of pt with eigenvalue ui a . The convenience 
of this choice arises in the time independence of the eigen¬ 
values u! a , which can be illustrated by the following ar¬ 
gument: 


As p t evolves according to Eq. (141, we must have 


d t Tr{p?} = Tv In 


^ [Hsys 




jn— 1 

Pt 


= 0 (16) 


for all positive integers n. Since the coefficients of the 
characteristic polynomial of pt can be expressed in terms 
of combinations of traces of powers of pt , it follows that 
the spectrum is time independent; i.e.; oj a {t) = w Q ( 0 ) = 


\Li ,01 ddsys] — -([ Lu],on Hsys]^ ~ w ^iL,a ( 12 ) It follows immediately that 


and so 

[Li, a Lui,p j H sys (\m = Li a [L u ,p,H sys (xm 

S' [Li, a i Hsys{X(t))]L Uj p 

= ljJ L l 

= 0 (13) 


establishing our claim. 

We see from Eq. (10 1 that the average information en¬ 
tropy is proportional to g 2 . This seems reasonable since 
if g = 0 , then there would be no coupling between the 


dtPt =^2^ad t [\uJa(t))(u}a{t)\] (17) 

a 

by the Spectral Theorem. Using the eigenket basis 
|w a (f)) affords us a simple expression for dtpt■ we need 
only compute the time-derivative of the projection op¬ 
erators since the eigenvalues w a are time-independent. 
However, the (time-dependent) energy eigenkets are more 
convenient for practical applications and so we express all 
relevant quantities in terms of this basis. 

Using the derivative truncation approximation, we can 
deduce the approximate eigenkets of pt in terms of the 
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energy eigenkets: 

K> 


“) +5Z 

b^a 


= l e a) + (<5|w a )) 


e -^a — e ~Ptb 

dX a 


d\ot 

(e b \S P M\e b ) — 


(18) 


A “ dt ' 

. This follows immediate from 


Furthermore, u a ~ - 
the approximation Eq. (15) and the definition of eigen¬ 
kets. 

The construction of the approximate eigenkets and 
eigenvalues forces the diagonal entries of Sp\a in the en¬ 
ergy eigenket basis to vanish. This is consistent with the 
requirement of positivity as both p t and p1 q are positive 
in Eq. (15), but the term involving 5p\<* is sensitive to 
the rate of change of A Q which could be negative. 

Using Eqs. © and , we see that 


dtpt K E Zt a\aP e “>M dt 


d\° 


(19) 


and so we obtain the expansion of Sp X a in the energy 
eigenket operator basis: 


a^c b 

hi 


/ hie 0eb 

V e ac Zt 


Tr{|e c )(e a |9A“ [|e 6 )(e 6 |]} 


= ^ J — ^p e t q ) ac \e a ){e c 


M(e c 

( 20 ) 


a^c 


with e ac = e a ~e c and (d x «p er! ) ac = Tr {\e c )(e a \d x <*Pt q }- 
Furthermore, from Eq. we have 

(*K>)». = £ e -„. Z l e -s.. ^z ( 21 > 

b^a 

By constructing an explicit expression for 8p x <* we 
may now construct a quadratic functional approximat¬ 
ing (E/) a : 

(S/) A = ff 2 / dt Tr{£(p t ) (-lnp t +lnp® 9 )} 

Jo 

C T 

dt —Tr{£(Sp xf ,)(-\np t +lnp e t q )} 

( 22 ) 

which follows from £(pt q ) = 0. 

We compute the trace in the |w a ) basis in which — In p t 
is diagonal. Since Tr{£ (•)} = 0, we have 

^-Tr {C(6p xf >) (-In p t )} 
dX 0 

= ^ Tr (<W) ln Pt + ln Z t)} 

~ ^-Tr{£(^ A /j)(-lnpj 9 )} 


dt 


+ Y / 2 ^aVl{(e a \C(5p x ,)6(\u a )) xa } 


dX a dX 0 
dt dt 


(23) 


From this follows a useful expression for the entropy 
production 


(^i)A^g 2 r 

Jo 


, dX a dX 0 

dt — -— x 

dt dt 

^213e a m{(e a \C{5p X f<)S (|w a )) AQ 

a 

+ {e a \C(5p x <*)8 (|w a )) A/3 } (24) 


after symmetrization in the a and f3 indices. We may 
write this expression as 


<£ 


I A 


dt 


E~, 

jklm 


dX a dX 0 


dt 


Ajkim {d x «Pt q ) jk ( d X ftp eq ) 
(25) 

where the components of the matrix A depend on C and 
the energy eigenvalues and eigenkets. We include the 
explicit components of A in the appendix. 


Eq. (24) approximates the average information entropy 


produced during a finite-time driving protocol of a quan¬ 
tum system weakly coupled to a large thermal bath us¬ 
ing only quantities directly calculable from the time- 
dependent system Hamiltonian and the Lindblad oper¬ 
ators. 


We will now use Eq. (24) to explicitly compute the av¬ 


erage information entropy produced by driving a simple 
two-state model quantum system. 


III. TWO-STATE MODEL SYSTEM 


We apply Eq. (24) to a simple two-state system with 
astern Hamiltonian 


H sys = - hw y (t)<Ty + hu z (t)cr z 


(26) 


where a y and a z are the Pauli spin matrices. Optimiza¬ 
tion of such a system could potentially be useful in ap¬ 
plications such as quantum annealing [47( I48J . 

We further assume that the system is coupled to a ther¬ 
mal bath of bosonic oscillators so that the full Hamilto¬ 
nian is 


Htot — H sys (t ) + g a z ® B -f Hb (27) 

where H B = u k b\b k , B = J2 k (b{ + b^j and g is a 
small coupling constant. 

It is convenient to work in an eigenbasis of a z repre¬ 
sented by the column vectors 


It} = 




In this basis, 


H S ys — ^ 


W z 

— lUJy 



(28) 


(29) 


lm 
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It is most convenient to express the control parameters 
in polar form: /3hco z = r cos(9), fihujy = rsm(O). In this 
form, it is not hard to show that the eigenvalues of H sys 
are r with eigenvectors 


Therefore, the optimal average entropy is given by 

<E ' U - “ ^ (.I ’ ° 



( — i sin (0/2) \ , . ( i cos (0/2) \ 

{ cos (0/2) { sin (0/2) )' 


(30) 


For simplicity, we assume r > 0 and 0 < 0 < f. A 
straightforward calculation demonstrates that 


r _/ COS 2 (60 * cos(0) sin(0) 

0 ' —icos(0) sin(d) — cos 2 (0) 


(31) 


|sin 2 (6 | ) — i cos 2 (0/2) sin(0) 

_2r — * sin 2 (0/2) sin(0) — lsin 2 (0) 


(32) 


^2r ~ 


lsin 2 (0) isin 2 (0/2) sin(0) 
„2ia /o'! l 


*cos 2 (d/2) sin(0) — |sin 2 (0) 


(33) 


are the operators defining C for this model system. 

Choosing the bath oscillator frequencies uik so that 
7(0) = 70 and 7(±2r) = 7ie ±r for 70,71 > 0 and setting 
h = 1, we have from Eq. (f2lj) 


If we choose Oq = 0 , Of = it /2, then 


(S/)y 


A 


70 (x 0 - \) In (^) 


2r 2 


EllipticE 1 — 


7l { e r + e~ r ) 
870 


(39) 


Eq. (39) exhibits the expected 1/r behavior of the op¬ 
timal average entropy production. Fig. [T] illustrates the 
dependence of (£/)A opt on the constants 70 and 71. We 
see that (T,j)\ opt has a 71 = constant profile described 
by s/lo- The quantities 70 and 71 are related to the 
bosonic frequencies uik of the thermal bath m and are 
thereby related to the noise of the quantum system. The 
bosonic frequencies consequently have a relatively sim¬ 
ple influence on the overall average entropy production 
in the finite-time long duration limit via Eq. (39). 


IV. CONCLUSION 


<£/>a 



7o 


(xo-l)ln(x^) 
2 r 2 


cos 2 (0) 


+ 71 



16r 2 


sin 2 ( 0) 





(34) 


where x 0 = : , e _ 2r . 

While this functional is non-negative as expected from 
the fluctuation theorem of pj], we see that it vanishes 
if 0 is held constant. If r is allowed to vary while 0 is 
fixed, no information entropy is generated on average. 
In terms of Riemannian geometry, this means that the 
metric tensor possesses null directions. 

We apply the Euler-Lagrange equation to obtain the 
optimal entropy production: 


implies 

(36) 

Upon integration, we find 

C 1 f 8f 

0=- C=- a[0}d0. (37) 

T Jeo 


Using the formalism of [TTTTl - ITTTTi . developed for classical 
systems, we were able to construct a general approxi¬ 
mation of the average information entropy of a quantum 
system driven in finite-time in terms of a quadratic func¬ 
tional of velocities in parameter space. This functional 
can be interpreted as endowing parameter space with a 
semi-definite metric in which optimal protocols are equiv¬ 
alent to geodesics. For a simple two-state driven quan¬ 
tum system weakly coupled to a thermal bosonic bath, we 
were able to derive an approximate expression for the av¬ 
erage information entropy. This expression has the char¬ 
acteristic 1/r dependence with a coefficient compactly 
expressed in terms of quantities related to the bosonic 
bath oscillator frequencies. 

Interestingly, in the simple two-state example studied 
here, the quadratic functional approximation for (£j)a 
we derived possesses null directions when expressed in 
terms of metric geometry on the space of control param¬ 
eters. In fact, for this model system, changing r results in 
a shift of energies. If r alone is changed, then it turns out 
that the system density matrix evolves exactly without 
a change in average information entropy. Consequently, 
our solution contains a null direction, which makes the 
quantum tensor not positive definite but positive semi- 
definite. It seems likely that this behavior is generic, 
though a proof is lacking. 

It is gratifying that this framework allowed us to ob¬ 
tain a general expression for the approximate entropy 
production in driven non-equilibrium quantum systems, 
as well as a closed-form solution for the minimum entropy 
production possible for a specific system with relevance 
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to quantum annealing. We are encouraged by the suc¬ 
cess of this approach for this simple system and we hope 
that this program will lead to further insight into the 
optimization of quantum systems out of equilibrium. 



FIG. 1. Plot of (Sj}A opt /(S/)A opt ,o where (S/)A opt ,o is the 
optimal average entropy for 70 = 1 , 71 = 0 . 
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VI. APPENDIX 


We record here the components of A appearing in 
Eq. (25) for the convenience of the reader. 


1 npz t 

2 “ (e~P eb - e~P e «) e dc 

b>a,d>c 


X { 


(fijd &kc $lb $ ma + Sm s me 3jb 4a) (e a \jC(\e d )(e c \)\e b ) 

“b (Sjc 4d $ lb & ma “b &lc $md &jb Ska ) (e a |£(|e c )(e d |)|e b ) 

T (Sjd Skc Sl a Smb -(- Sid S mc Sj a Skb ) {^b |£( | ^d) (^c |) |^a) 

(Sjc Skd Sl a Smb ~\~ Sic Smd Sj a Skb) {^b |£( | ^c) {^d |) |^a) } 

(40) 
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